Quasiparticle poisoning and Josephson current fluctuations induced by Kondo 

impurities 
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We introduce a toy model that allows us to study the physical properties of a spin impurity 
coupled to the electrons in the superconducting island. We show that when the coupling of the 
spin is of the order of the superconducting gap A two almost degenerate subgap states are formed. 
By computing the Berry phase that is associated with the superconducting phase rotations in this 
model, we prove that these subgap states are characterized by a different charge and demonstrate 
that the switching between these states has the same effect as quasiparticle poisoning (unpoisoning) 
of the island. We also show that an impurity coupled to both the island and the lead generates 
Josepshon current fluctuations. 
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Introduction.— Superconducting circuits based on 
small Josephson junctions are promising candidates for 
the implementation of qubits [l|, 0, d, U, [B| and for the 
development of a prototype quantum current standard 
Unfortunately, the performances of these devices 
arc significantly limited by different types of noise whose 
sources remain mostly unknown. Particularly danger- 
ous one for Single-Cooper-pair Transistors and Cooper 
pair boxes is the noise produced by the incoherent tun- 
neling of single quasiparticles into the superconducting 
island. Each tunneling event changes the island charge, 
thereby shifting the operation points of the device. An 
important requirement for the regular operation of these 
devices is that this tunneling is very rare. Despite signif- 
icant experimental efforts to reduce quasiparticle poison- 
ing @, H, H, 03, 11 , 12] a complete understanding of its 
microscopic mechanisms is still missing. The goal of this 
Letter is to show that the mechanism of the charge noise 
discussed in Ref. [l3| might be also responsible for the 
creation of the low energy quasiparticle traps and pro- 
vides an explanation of the puzzling features observed in 
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The work [13[ shows that Kondo-like traps located at 
the Superconductor Insulator (SI) interface might pro- 
duce the charge noise in small Josephson charge qubits; 
similar mechanism might be responsible for the critical 
current fluctuations in large superconducting contacts 
[l5| . These Kondo traps arc impurities with a singly 
occupied electron level that carry a spin degree of free- 
dom. Each trap is characterized by an effective Kondo 
Temperature Tk that depends exponentially on its hy- 
bridization with the conducting electrons in the bulk su- 
perconductor. In this mechanism both charge and crit- 
ical current noise originate from the electrons tunneling 
between those Kondo traps with Tk ~ A, where A is 



the superconducting gap. In this Letter we show that 
Kondo-like traps might also be responsible for quasipar- 
ticle poisoning of the superconducting island. Further, 
we show that such traps located close to the Josephson 
junction generate additional sources of critical current 
fluctuations due to their coupling to the superconductors 
on both sides of the barrier. This mechanism for critical 
current noise provides the alternative to the conventional 
picture of fluctuators blocking conducting channels in the 
insulating barrier. In order to derive these results we in- 
troduce a toy model that captures the essential physics 
of a spin impurity coupled to the superconducting elec- 
trons in the superconducting island. By computing the 
Berry phase that is associated with the superconducting 
phase rotations in this model, we show that two differ- 
ent low energy states of the impurity are characterized 
by a different charge. As a consequence, switching be- 
tween these two low energy states has the same effect as 
quasiparticle unpoisoning (poisoning) of the island. Fi- 
nally, we use this model to study the effect of the motion 
of electrons between the Kondo-like traps in a Josephson 
junction and we prove that if one of those traps is coupled 
to both the lead and the island, these processes result in 
critical current fluctuations. We begin with the review 
of the features of the Kondo physics that are relevant for 
the following and which provides justification of the toy 
model. 

The behavior of a spin- 1/2 impurity coupled antiferro- 
magnetically with an exchange constant J to an electron 
gas characterized by a constant density of states po within 
a bandwidth D is completely different at high and low 
temperature regimes. In the former, the electrons scat- 
ter off the impurity inelastically in a spin-flip process 
while at low temperatures the impurity is screened by 
the electrons forming a bound singlet state leaving only 
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clastic scattering. The crossover takes place at energy 
scale of Kondo temperature Tk ~ De~ 1 ' Jp °. All rele- 
vant physics is described by the Anderson Hamiltonian: 

lead 



H = Hi Pnd + H d + H sd , where Hi ead = e K c* c K a and 



H d = edc] la c d a + Uridtndi 

a 

H s d = ^ ( Vc i<r c d(7 + h.c.) 



(1) 



Here c Ka {c\ a ) is the annihilation (creation) operator for 
an electron in the band, Cda (c d <j) ^ s ^ ne annihilation (cre- 
ation) operator for an electron on the impurity site. U 
denotes the Coulomb onsitc repulsion of the trap and 
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Cda- In the limit of strong onsitc repulsion, i.e. 
U r, where T = irpoV 2 defines the hybridization to the 
conducting electrons, using the Schrieffer- Wolff transfor- 
mation [16( , we can map the Hamiltonian given in Eq. {]]) 
to the Kondo model [171 ]: 



H = H, 
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(2) 



Here J KK > = J = 8V 2 /U. This mapping neglects the par- 
ticle hole asymmetry of the original problem, the effects 
that are small in Tk/D but might have important phys- 
ical consequences. Notice that in the limit J/D — > 1, the 
Kondo temperature becomes Tk ~ J- 

A more complicated problem is presented by the im- 
purity interacting with the conduction electrons in the 
superconductor. At present, only qualitative arguments 
and numerical results are available which show that 
the competition between the Kondo temperature of the 
trap and the superconducting gap A of the lead re- 
sults in three different regimes for the system (impu- 
rity+superconductor): (i) Tk <C A, where the ground 
state of the system is a doublet and it is characterized 
by an odd number of electrons; (ii) Tk 3> A where the 
ground state of the system is a singlet, the electron of the 
impurity forms a bound state with the superconducting 
electrons and the total number of electrons is even; (iii) 
Tk ~ A where the singlet and doublet states become al- 
most degenerate. 

Toy model and quasiparticle poisoning.— In order to 
formulate a simplified model that captures the effects of 
an impurity interacting with the superconducting elec- 
trons, we notice that the Kondo physics can be viewed 
as a result of the 'poor-man scaling' in which the high 
energy degrees of freedom are gradually integrated out 
resulting in the logarithmic growth of the effective inter- 
action. In the presence of the superconducting gap the 
process of integration has to stop at A. If at this moment 
the interaction is comparable with A, a bound subgap 
state can be formed in agreement with the results de- 
scribed above. This shows that the essential physics can 
be captured by a simple model in which the spin interacts 
with a single electron mode with the coupling constant 



J ~ A and energy e < A that is described by the BCS 
Hamiltonian: 
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The interaction between the impurity and the supercon- 
ducting electrons is described by the spin exchange cou- 
pling given in Eq.([2|): 

Htoy = Hbcs + Jk.k'S ■ cj co . i fT (Tl(T2 c K ' (J2 (3) 

k,k' 

We assume that the coupling is isotropic: i.e. 
Jk.k' = J ~ Tk ~ A. The Hamiltonian ([3]) can be 
readily diagonalized. We choose the state basis: 
{\i)kh \j)-ki, \<y)i mp } where i,j = 0, 1 denotes respec- 
tively the absence or presence of the quasiparticle in the 
single electron mode while a =f|",JJ- represents the spin 
configuration up or down of the electron in the trap and 
find the lowest eigenvalues: 



E Q = £ - v/A2 

E 1 =e- \t k 



(4) 



corresponding respectively to the (non normalized) dou- 
blet and singlet states: 
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As expected, the spin impurity interacting with the con- 
duction electrons in the superconductor leads to the for- 
mation of weak Kondo subgap states. Notice that the 
subgap states have different properties: the doublet is 
characterized by an odd number of electrons, its degen- 
eracy is due to the spin degree of freedom of the trap 
while the singlet state is a maximally entangled state 
with even number of electrons. Depending on the ra- 
tio Tr-/A, the ground state of the system can be either 

doublet or singlet. At a special value T K = — \J e 1 + A 2 

singlet and doublet states are degenerate while for traps 
with Tk ~ T K ~ A singlet and doublet states are almost 
degenerate. 

So far, we have reproduced in a simplified manner the 
main content of the numerical results. We want now 
to show that singlet and doublet states differ in another 
important feature: namely, the doublet corresponds to 
a non zero off-set charge in the superconducting island. 
For this purpose, wc recall that the operator n, describing 
the excess number of Cooper pairs on the island and the 
operator 9, representing the superconducting phase, are 
conjugate variables, i.e. [0,n] = i, and that the Hamilto- 
nian of the island/box H is i/ box is invariant with respect 
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to the local gauge transformation: U 1 H is i/i )ox U where 
U = e" LgB and n g is the off-set charge induced in the is- 
land in unit 2e, that plays a role similar to the vector 
potential appearing in the Hamiltonian of an electron in 
a magnetic field. One consequence of these observations, 
is that we can deduce the value of the off-set charge in- 
duced by the Kondo impurity in the singlet and doublet 
states from the value of the Berry phases associated to 
these states. We find that 

2 71 
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Notice that a Kondo impurity induces a non-zero off-set 
charge in the superconducting island only when the sys- 
tem (impurity+superconductor) is in the doublet state, 
while in the singlet state the electron charge is fully 
screened. Moreover, when e — ► 0, we find that n — > ^, i.e. 
exactly one electron is induced on the superconducting is- 
land. Let us now consider two Kondo-traps with Tk ~ A 
located at the SI interfaces, one on the superconducting 
island, another on the lead within distance £ from each 
other. The electron tunneling process across the junc- 
tion couples these traps. When the Kondo trap located 
at the SI interface in the island switches between singlet 
and doublet, the parity of the island changes from even 
to odd. For this process to be physically relevant, the 
energy difference between these states should be smaller 
than T. Thus, the pairs of Kondo subgap states with 
close energy levels might be responsible for the quasipar- 
ticle poisoning in superconducting devices. Notice that if 
the two traps are located on the same side of the barrier, 
the switching between singlet and doublet caused by the 
tunneling of quasiparticles through the superconductor 
results in charge fluctuations, i.e. the entire process can 
be viewed as a new type of a charge fluctuator. 

We now discuss the implications for the recent ex- 
periments where quasiparticle tunneling rates were mea- 
sured with microsecond resolution [lCJ, 
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micrometer-sized island with capacitive gate electrode 
that was probed by two Josephson junctions. The 
island charging energy is modulated by the gate as 
E?(n g ) = E c (n-n g ) 2 i where E c = e 2 /2C, e is the 
electron charge and C is the total island capacitance, 
n g = C g Vg/e is the normalized gate charge and n is the 
integer number of excess charges on the island. At n g = 1 
the electrostatic energy of the system is minimized when 
unpaired electrons reside in the superconducting island. 
Thus, at n g = 1 the island is a trap for a quasiparticle 
with depth SE = E C - Ej/2 + A; - A*. Here A ?; , A; are 
the superconducting gap of the island, lead and Ej is 
the Josephson energy. A model suggested by Aumen- 
tado et al. Q explains many features of quasiparticle 
poisoning of the island. In this model some unknown 



non equilibrium source of quasiparticles produces them 
in the leads. Quasiparticles are able to tunnel onto the 
island which acts as a trap. Subsequently the trapped 
quasiparticle is thermally excited (unpoisoning) out of 
the trap and the island returns to its even state. Rele- 
vant implications of this model are that the quasiparticle 
poisoning can be reduced by putting normal metal leads 
(QP traps) close to the junctions in order to filter the 
quasiparticles and by making A, > A;, because it works 
as a barrier, which prevents non equilibrium quasiparti- 
cles in the leads from entering the island. Experiments 
showed that these ideas help to reduce quasiparticle poi- 
soning, but do not eliminate it. In particular, the effect 
of quasiparticle traps have been recently studied in Ref. 
[H |. In these experiments, two similar Cooper pair boxes 
were fabricated with (QT) or without (NT) quasiparti- 
cles traps attached to the leads. The island was biased at 
rig = 1 and the dynamics of the quasiparticles captured 
by the island was characterized by their incoming (t even ) 
and outgoing (t 0( id) rates. One expects that the incoming 
process involves quasiparticle tunneling into the island 
and its relaxation to the bottom of the well, while the re- 
verse process involves thermal excitation. As a result, the 
outgoing rate should be smaller by a factor cx e - SE / KT 
than the incoming rate. This is in contrast with the data 
that show that the trapping and escape rates are roughly 
equal and temperature independent below T < 200mK. 
However, their values are dramatically different in the de- 
vices with or without traps: t even ~ t a dd ~ (10 2 — 10 3 )/is 
(QT) and t even ~ t odd « (0.1 - l)/is (NT). 




FIG. 1: Sketches of quasiparticle poisoning and unpoisoning 
due to quasiparticle tunneling between weak Kondo subgap 
states located at the lead/island SI interfaces. 

The presence of subgap states in the lead and the is- 
land provides a different scenario where two new pro- 
cesses are present: quasiparticles with energies above 
Aj — E c + Ej/2 tunnel from a subgap state in the lead 
to the continuum in the island while quasiparticles below 
these energies tunnel between the subgap states in the is- 
land and in the lead (see Fig. [I}. The rate of the former 
process is T sc = GS, where G ~ 1 is the conductance of 
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the barrier in the units of e 2 /h and S = 1/ViVAi is the 
typical level spacing in the island. For a typical island 
of volume Vj = 750nm x 125nm x 7nm and a typical Al 
electron density of states vai ~ 35/eVnm 3 , we estimate 
r sc w 10 _7 s _1 . The rate of the exchange process between 
subgap states is much slower T ss <C r sc because it occurs 
between two localized states and it depends on the level 
width of the state. In both cases the tunneling does not 
involve a significant energy transfer, so we expect these 
rates to be temperature independent. In the presence of 
quasiparticlcs with energy larger than Aj — E c + Ej/2 
the first process dominates and the observable rate is 
T sc . The presence of quasiparticle traps attached to the 
leads eliminate high energy quasiparticles and the rate 
decreases to T ss in agreement with observations. This 
scenario can be checked by fabricating QT devices with 
traps in the lead located at different distance L ~ £. 
The presence of these traps at sufficiently small distance 
broadens the levels of the subgap states in the leads and 
consequently it should increase the tunneling rate into 
the island as cx e~ L ^, where £ is the coherence length 
of the superconductor. The estimates of the rate assume 
that only a few subgap states are active at the same time, 
large number of these states would make the effective rate 
higher; we do not know the density of these states and 
their occupation in a realistic system. 

Josephson current fluctuations.— Kondo-like traps 
with Tk ~ A provide additional source of the critical 
current fluctuations when they are located close to the 
Josephson junction barrier because in this case the spin of 
the Kondo trap is coupled both to the electrons in the is- 
land and in the lead. This can be easily seen by including 
in our toy model the electrons of the lead and a tunnel- 
ing between the lead and the island. The Hamiltonian 



becomes H = H 



(i) 

toy 



H 



(2) 
BCS 



Hrj? and the quasiparticle 



tunneling through the junction barrier is given by: 



T, 



4it C *aT + C - 



c-k 2 i + h.c. 



We assume that the superconducting leads are equal and 
we calculate the correction to the lowest energy eigenval- 
ues and eigenvectors at the second order in perturbation 
theory in the tunneling |Tfe li fe 2 | ~ T. We find that the 
correction to the singlet state depends on the phase dif- 
ference ip = 9\ — #2- The dependence on <p implies addi- 
tional contribution to the Josephson current. A straight- 
forward but lengthy calculation, gives the contribution 
of the Kondo impurity in the barrier to the Josephson 
current: 



Sir. 



J e 2 



■ Sill if 



(6) 



To find the parameters of the toy model that describe 
the physical situation in which a many channel junction 
couples a small metallic island to a large superconducting 
lead, we compare the pairing field induced on the state 



ki in the islands by the superconducting order in the 
leads in the realistic situation {yT 2 ~ GS) with the cor- 
responding quantity in the simplified model (T 2 /A) and 

get — = GS. Notice that this reasoning holds for small 

superconducting islands whose size is less than the super- 
conducting correlation length £. For larger islands the 
impurities at a distance larger than £ from the junction 
are coupled exponentially weakly to the superconductor 
on the other side of the barrier. 

Conclusions.— We have shown that subgap states gen- 
erated by magnetic impurities due to the competition 
between superconducting pairing and Kondo effect act 
as very efficient quasiparticle traps. We argued that 
the presence of such states in a typical Single-Cooper- 
pair Transistor and Cooper pair box might explain the 
results of recent experiments where unexpected poison- 
ing/unpoisoning rates were observed. We have also 
shown the same subgap states generate critical current 
noise. 
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tract number W911NF-06-1-0208 and NSF ECS 0608842. 
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